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Abstract 

In one dimension the coupling of electrons to phonons leads to a transition from a metallic to 
a Peierls distorted insulated state if the coupling exceeds a critical value. On the other hand, in 
two dimensions the electron-phonon interaction may also lead to the formation of Cooper pairs. 
This competition of superconductivity and charge order (in conjunction with a lattice distortion) is 
studied in this letter by means of the projector-based renormalization method (PRM). Increasing 
the electron-phonon interaction, we find a crossover behavior between a purely superconducting 
state and a charge-density wave where a well-defined parameter range of coexistence of supercon- 
ductivity and lattice distortion exists. 
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Introduction. In recent years, the interest in two-dimensional electron-phonon models has 
been considerably renewed, mainly triggered by strong experimental indications that the 
electron-phonon interaction has a substantial impact on the properties of high-T c cuprates 
and by the discovery of superconductivity in MgB 2 below a rather high T c of about 
39 K [2]. Furthermore, in the context of quantum-phase transitions, such systems are also 
of general theoretical interest because the electron-phonon interaction causes a competition 
jetween a structural instability and superconductivity. As pointed out by many authors (e.g. 
2I, Q], 0, 0]), this competition may play an important role in Alb materials, (Ba,K)Bi0 3 
or high-temperature superconductors. 

The simplest general model of a coupled electron-phonon system on a square lattice is 
given by 

U = Ho + Hx (1) 

k,cr q 
Hl = 9* { 6 q C L C k+q,- + &qC k+qi(T C ki(T } . 



k,q,cr 

Here, , Ck, CT and &k are the creation and annihilation operators of electrons and 
phonons. Assuming an electron hopping between nearest-neighbor sites on a square lat- 
tice, the electronic dispersion is given by = — 2t(cos k x a + cos k y a) — //, where [i is the 
chemical potential. Moreover, cu q is the phonon energy, and g q is the wave vector dependent 
coupling between the electrons and phonons. In the case of dispersion-less phonons c^o and a 
q-independent coupling g, the model ([1]) reduces to the Holstein model. Then, the last term 
describes a local interaction between a lattice displacement and the local electronic density. 
For a more realistic description of high-T c materials, we would have to add a Hubbard-U 
term to the model. 

The interaction between a possible superconducting (SC) state and an insulating Peierls 
charge- density wave (CDW) phase in the 2d-Holstein model has been subject to a num- 
ber of Monte Carlo (QMC) studies [zj|. However, their mutual influence is not yet fully 
understood, because only correlation functions, namely SC and CDW susceptibilities, can 
be studied within QMC calculations but not the fundamental order parameters themselves. 
Nevertheless, it is clear that the two phases compete: the appearance of strong correlations 
of one kind suppresses the development of correlations of the other kind. But these QMC 



studies could not clarify whether a coexistence of the two phases is possible or not. 
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On the other hand, particular attention received a work by Bilbro and McMillan 
1976 where the competition of superconductivity and the martensitic transformation in A15 
compounds was studied by a mean-field treatment of an appropriate model Hamiltonian 
and where a possible coexistence of a Peierls gap and a SC gap for a common portion of 
the Fermi surface was proposed. The Bilbro-McMillan model has been the basis for further 
theoretical work till now, for instance in the context of high-temperature superconductors 

in n 

9|],[10( and heavy fermion compounds [11]. However, to our knowledge two open problems 
in the Bilbro-McMillan model remain to be solved: First, all fluctuations beyond mean- field 
theory are neglected so that the equations for SC and CDW gaps are given by weak coupling 
expressions and any influence of a possible renormalization of the phonon energies is not 
included. Second, the electron-phonon interaction should be part of the basic microscopic 
model to have access to phonon properties. 

Here, we tackle these two open questions by means of the recently developed projector- 
based renormalization method (PRM) [3]. In this way we are able to include fluctuations 
beyond mean-field theory and to take into account additional renormalization effects like 
phonon softening. Furthermore, we directly access the order parameters of the two phases 
where, in agreement with reference [8(, we find a parameter range of coexistence of SC and 
CDW at half-filling. Depending on the strength of the electron-phonon coupling the Peierls 
phase will be suppressed by the superconducting phase or vice versa. 

Projector-based renormalization method. The PRM 12] starts from the usual decom- 
position of a many-particle Hamiltonian into a solvable unperturbed part Tio and a per- 
turbation Tii where Tii does not contain any part that commutes with 7{q. Thus, the 
perturbation Ji\ consists of transitions between the eigenstates of 7io with non-vanishing 
transition energies. The basic idea of the PRM is the construction of an effective Hamil- 
tonian 7i\ = 7^o,a + T~ti,\ with renormalized parts ?io,\ and TCi } \ where all transitions with 
transition energies \Eq x — E™ x \ larger than a given cutoff energy A are eliminated. Here, 
Eq X and E™ x denote the eigenenergies of 7^o,a- 

The renormalization procedure starts from the cutoff energy A = A of the original model 
7i and proceeds in steps of AA to lower values of A. Every renormalization step is performed 
by means of an unitary transformation where the elimination between the cutoffs A and 
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(A — AA) reads 



W(a-aa) = e ' AA 7~C\ e~ ,AA . (2) 

Here, the generator Xa,aa of the unitary transformation has to be fixed appropriately (for de- 
tails see Ref. Il2l). In this way difference equations are derived which connect the parameters 
of H\ with those of TC(\-a\), an d which we call renormalization equations. 

The limit A — > provides the desired effective Hamiltonian H = Hx^o = 7^q,a->o where 
the elimination of the transitions originating from the perturbation Tlx leads to a renormal- 
ization of the parameters of 7^o,a^o- Note that 7i is diagonal or at least quasi-diagonal and 
allows to evaluate physical quantities. The final results depend on the parameter values of 
the original Hamiltonian TC. Finally, note that Ti. and 7i have the same eigenvalue problem 
because both Hamiltonians are connected by an unitary transformation. 

SC and CDW phases in the half-filled Holstein model We now want to apply the PRM 
approach to the Holstein model at half-filling in order to study the interplay between SC 
and CDW phases. For this purpose a uniform description of both the SC and the insulating 
CDW phase has to be found. In the SC state the gauge symmetry is broken. Therefore, 



following Ref. [13j, a field which breaks gauge symmetry should be added to the Hamiltonian. 
Similarly, for half-filling one has to take into account that the unit cell of the system can be 
doubled in the case of a dimerized insulating ground state jl4j ]. Therefore, we add to the 
Hamiltonian symmetry breaking fields as follows 

n H + E ( A k c L Aa + Afc-k^) (3) 

k 

+\ E ( A k c k,<T C k-Q,o- + h.c) + VNA b (bl + b Q ) 

k,cr 

where all fields are assumed to be infinitesimally small (A^ — > 0, A£ — > 0, A b — > 0). Q = 
(n/a, it /a) is the characteristic wave vector of the CDW phase. Due to the doubling of the 
unit cell in the insulating phase the Hamiltonian is best rewritten in the reduced Brillouin 
zone where we have £k-Q = — £k- Exploiting in addition the coupling of the creation operator 
c ic a t° * ne annihilation operator c_k.-o- due to superconductivity a four- dimensional compact 
vector notation for the electronic one-particle operators can be introduced 

= ( C-k+Q,j c m c -k,l c Lq,t ) ( 4 ) 
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The renormalized Hamiltonian, after all transitions with energies larger than A have been 
integrated out, can again be divided into 7i\ = 7io y \ + Hi, a- For 7Y ,a we make the ansatz 



n , x = *£^£*a 



(5) 



ker.BZ 



+ E E 

o=0,l qGr.BZ 



with 



H 



k,A 



/ e -A p A s *\ 
e A s A p 

_ A P A S* _ £ Q 

\ A s A p -e / 

x / k,A 

where all parameters now depend on k and A due to renormalization processes. Note that 
the symmetry breaking fields have been included in 7io,A ; and the phonon field term ~ A 6 



(t) 

0,q 



M) 
q-Q 



from (j3J) was incorporated in redefined phonon operators. The new operators fa, 
and fej q = bq characterize the phonon branches in the reduced Brillouin zone with energies 
^o,q,A = ^ q -Q,A and ^i, q ,A = ^q,A- The renormalized interaction TCi ; \ can be compactly 
written as well 



Hl * = ~m Z Z) { 6 lq^^"k,q,AC k+q + h.C.} 



(6) 



k,qGr.BZ 7 e{0,l} 



where the elements of the new 4x4 matrix Hi 



1 k A again depend on A and on wave vectors 
q and k. Note that the general structure of TCq^ an d TCi y \ and thus of the matrices H^x an d 



H 



7,k,q,A 



remains always the same during the renormalization procedure and agrees with that 



of the corresponding matrix of the general electron-phonon model (pQ). The initial values of 
the A-dependent parameters in (jSJ), (jSJ) are fixed by the original Holstein model. 

The eigenvalue problem of 7^ ,a can be solved analytically. For this purpose we introduce 



new A dependent operators a a ,k,x, (a = 1. 
vector notation is used. Here, we defined 

/ u 


Ac.A 



A)'- o-kA = -^k,AC k where the four-dimensional 



_,,p 



v p -v s \ 




u —v 

W F V s u 



\ V s 



V 1 







(?) 



u J 



k,A 
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The condition (wk,A) 2 + (^a) 2 + (^a) 2 = 1 guarantees that the usual anti-commutator 
relations for a aj k,A are fulfilled, a) a k A , cv^'.a = ^k,k'^a,a'- For the electronic part of Hq,a 



one obtains 



^A - ^ \ El ^ ( a l,k,A a l,k,A + 4,k,A a 2,k,A) 
ker.BZ 

+^2,k,A (^ a 3,k,A a 3,k,A + a 4,k,A a 4,k,A j j (8) 

where the energies are given by 



£iAM = ±V4A + l A k,Al 2 + |Ak,Al 2 (9) 

fc> r ^k,A > 0, whereas for e^x < the ±-signs have to be reversed. 

In order to find the renormalization equations (which governs the A dependence of the 
parameters of 1~C\), one has to evaluate the unitary transformation (j2j) explicitly. We use 
the following ansatz for the generator X\ t A\ 

X XA x = (10) 

= 4^ 2 {H,q4,A^£ q A ak +q ,A-h.C.} 
ViV k,qGr.BZ 7 G{0,l} 

Besides the fact that the generator has to be anti-hermitian, X\ AA = — Xa,aa, the operator 
structure of the ansatz (Tl0|) agrees with that of the interaction TC^x of ([6]) where the electronic 
operators c k are replaced by the eigenmodes a k A of 7io,A- The matrix v4^' k q in fflOl) has to 
be fixed in such a way that, with respect to 7^o,(a-aa)> only excitations with energies smaller 
than (A — A A) contribute to TCi^x-ax)- The renormalization equations for u; 7iqi A and the 
parameters of H^ x an d H^ k <i X are obtained by comparing Eqs. (J5j) and ([6]) with the result 
of the explicit evaluation of (T5]) using ansatz ( flOl) for Xa,aa, after all fermionic creation 
and annihilation operators a k have been transformed back to the original operators c k . To 
evaluate Eq. ([2]) an additional factorization approximation must be employed in order to keep 
only operators of the structure of those of ([5]) and ([6]). Therefore, the final renormalization 
equations still depend on unknown expectation values. They are best evaluated with the 
full Hamiltonian 7i in order to incorporate important interaction effects [l^j]. Therefore, 
we have to apply the sequence ([2]) of unitary transformations also to operators, -4(a-aa) = 
e Xx < AX Axe~ Xx,AX and exploit (A) = (Ax)n x - This procedure is performed for the fermionic 
and bosonic one-particle operators, cjp and by)q, where the same approximations are used as 



for the Hamiltonian. The final set of coupled renormalization equations is solved numerically. 
Thereby, the equations for the expectation values are taken into account in a self-consistency 
loop. The fully renormalized Hamiltonian is obtained for A — > 0, where the interaction Tii 
is completely used up for the renormalization of the parameters of 7io,A^o = W. Thus, an 
effectively free model is obtained, 

W = E + E E ^,q&U&r, q + £, 

ker.BZ qer.BZ 7 g{0,l} 

(11) 

where i^ 1 = H^ x _^ has the form as in Eq. (JSj). The final Hamiltonian 7i is diagonal 
and can be used to investigate various physical properties. Note, in particular that the 
eigenenergies £a/2,k and u) 7iq of 7i can be considered as quasiparticles of the full Hamiltonian. 

Two aspects of the presented PRM approach should be noticed at this point: Even 
though a factorization approximation must be employed in order to derive the effectively 
free model (TTTT) . fluctuations beyond mean-field theory are taken into account due to the 
renormalization procedure. Furthermore, both electron-phonon interaction and phononic 
degrees of freedom are included in the microscopic model so that a direct access to phonon 
properties is provided. 

Results and discussion. We consider a square lattice with N = 144 sites and restrict 
ourselves to the Holstein model with dispersion-less phonons c<j q = ujq. A local electron- 
phonon coupling g q = g is assumed. The temperature is set equal to T = and a small value 
of lvo is chosen, ujq = O.lt. In the following we concentrate on the s-wave like superconducting 
pairing because we could not find any stable ci-wave solution. Panel (a) of Fig. [1] shows 
the order parameters for the Peierls state, P = 4p Xlk( c k o- c k+Q,o-) (black), and for the 
superconducting state, S = J2k( c li,i c -k,i) ( re d) as function of g, where g is restricted to 
small values g/2t < 0.04. For comparison, in panel (b) of Fig. [I] the k dependent symmetry 
breaking fields A£ (black) and A^. (red) are shown for k = (tt/2,tt/2) for the same parameters 
as in panel (a). Note that due to Qj, A£ and A^. contribute together to the energy gap in 
the quasiparticle spectrum of H and either P and S or A£ and A^ can be considered as 
order parameters. As can be seen, for small values of g/2t < 0.010 the system is in a pure 
superconducting state, i.e. no charge order is present (P = 0). The superconducting gap 
increases roughly proportional to g 2 . In the intermediate g range, 0.010 < g/2t < 0.023, a 
coexistence of both order parameter P and S is found. Thus, the system is in a combined 
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FIG. 1: (Color online) (a) Order parameters P = f Ek^L^+Q,") (black line) and S = 
jj X^k( c k t c -k i) ( re< ^ nne ) as function of the electron-phonon coupling g for a square lattice with 
144 lattice sites at half-filling, uo/t = 0.1 and T = 0. (b) Renormalized values of the Peierls gap 
A£ (black line) and of the superconducting gap Aj^ (red line) at wave vector k = (n/2, 7r/2) for 
the same parameter values as in the upper panel. 

superconducting-charge ordered state. The g dependence of A s is no longer quadratic as 
in the small g regime. Instead A s reaches a maximum value and drops down to zero with 
increasing g. For g/2t > 0.023 the superconducting phase is completely suppressed and the 
system is in a pure charge ordered state. 

The mutual influence of the two order parameters is considered in Fig. [2j First, in panel 
(a) of Fig. f2]the renormalized superconducting energy gap A s is shown as function of g for two 
cases: (i) the former result from Fig. [1] which follows from the full renormalization equations 
(in red), and (ii) the result when the charge order is suppressed in the renormalization 
equations 'by hand' (in black). The comparison shows that superconductivity becomes 
strongly suppressed when the charge order is present for large g values. 

In panel (b) of Fig. [21 the dimerized lattice displacement (&q + 6q) is shown as function 
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FIG. 2: (Color online) (a) Renormalized superconducting gap as function of g for the case that 
CDW order is suppressed by hand (black curve). The red curve shows the complete solution 
including CDW order, (b) Dimerized lattice distortion (6q + &q) as function of g. In the black line 
a possible superconducting phase is suppressed The red line shows the complete solution allowing 
also a superconducting state. 

of g for the same parameter values as in Fig. [TJ Note that a dimer-like shift of the ionic 
equilibrium positions comes always along with the presence of a charge density wave. Thus, 
the expectation values (6q + 6q) are nonzero in the Peierls phase and can be considered as 
an alternative order parameter. The red curve in panel (b) of Fig. [2] shows (bL + 6q) as 
it follows from the solution of the full renormalization equations. For comparison, in the 
black curve the superconducting phase is artificially suppressed. The comparison of the two 
curves shows that for weak coupling g/2t < 0.01, when the superconducting state exists, 
the Peierls state is strongly suppressed. However, for larger values of g the superconducting 
order parameter vanishes and the red curve converges to the black curve. For the case, when 
the superconductivity is suppressed, the ionic shift (6j, + 6 Q ) increases almost linearly with g 
and a Peierls state is found up to g — > 0. In this respect the two-dimensional Holstein model 
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differs from the one-dimensional Holstein model, where a Peierls state exists only above a 
critical electro n-phonon coupling g c > 0. 

Summary. The developed approach to the two-dimensional Holstein model has two par- 
ticular advantages: Fluctuations beyond mean- field theory are taken into account, and the 
appearance of superconductivity and CDW state is directly proven on the basis of their 
order parameters. Thus, the presented PRM approach overcomes important limitations of 
QMC studies {t| and of the Bilbro-McMillan model and our results provide a reliable 
proof that indeed superconductivity and lattice distortion coexist in a certain parameter 
range of the electron-phonon coupling. 
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